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THE SPECTRUM OF THE DIFFUSE GALACTIC LIGHT I: THE MILKY WAY IN SCATTERED LIGHT 

Timothy D. Brandt and B. T. Draine 1 

Draft version November 18, 2011 

ABSTRACT 

We measure the optical spectrum of the Diffuse Galactic Light (DGL) - the local Milky Way in 
reflection - using 92,000 blank-sky spectra from the Sloan Digital Sky Survey (SDSS). We correlate 
the SDSS optical intensity in regions of blank sky against 100 fim. intensity independently measured 
by the COsmic Background Explorer (COBE) and InfraRed Astronomy Satellite (IRAS) satellites, 
which provides a measure of the dust column density times the intensity of illuminating starlight. The 
spectrum of scattered light is very blue and shows a clear 4000 A break and broad Mg b absorption. 
This is consistent with scattered starlight, and the continuum of the DGL is well-reproduced by a 
simple radiative transfer model of the Galaxy. We also detect line emission in Ho, H/3, [N n], and [S n], 
consistent with scattered light from the local interstellar medium (ISM) . The strength of [N n] and 
[S n] , combined with upper limits on [O in] and [He i] , indicate a relatively soft ionizing spectrum. We 
find that our measurements of the DGL can constrain dust models, favoring a grain size distribution 
with relatively few large grains. We also estimate the fraction of high-latitude Ho which is scattered 
to be 19 ±4%. 

Subject headings: ISM: dust, extinction; Scattering; Methods: statistical 



1. INTRODUCTION 

All astronomical observations include light from diffuse 
sources other than the target object. In ground-based 
data, the dominant sources of contamination in the opti- 
cal are airglow, scattered sunlight, and artificial sources. 
Space-based missions must still contend with zodiacal 
light, and all observations will include some emission and 
scattering from the Galaxy's interstellar medium (ISM) 
- the diffuse Galactic light (DGL). 

The first quantitative measu rements of the DGL w ere 
photoelectric measurements bv lElvev fc RoacH (|1937f) at 
A ?» 4500 A; after subtracting the zodiacal and air- 
glow contributions, the DGL was detec ted for |b| < 35° . 
Subs equent studies from th e ground (lElsasser fc Hau 
1960D and a sounding rocket (fWolstencroft fc Roselll96 l 
found intensities at |6| ~ 5° corresponding to ~ 50 
10 th magnitude stars per square degree, or XI\ « 
5 x 10~ 4 ergc m~ 2 s~ 1 sr 1 . Sub s equent obseryations 
from satellite s dLillie fc Wittl 119761: [Morgan et all 1 197: 
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tended these studies into the vacuum ultr aviolet (UV). 
Most of these studies were broadband, but iMartin et all 
(1990) detected fluorescent emission from UV-pumped 
H 2 . 

The spectrum of the DGL contains a wealth of informa- 
tion about the physical environment where it originates 
and the dust that emits or scatters it into our line-of- 
sight. In this paper, we present a novel way of measuring 
the spectrum of light scattered by the Galactic ISM. This 
is a spectrum of the Galaxy in reflection, plus possible 
luminescence from interstellar dust. 

2. METHODOLOGY 
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The optical surface brightness of the DGL, XI \ as 

10~ 4 ergcm~ 2 s _1 sr , is far too low to measure a spec- 
trum directly. Further, any such spectrum would be 
a combination of terrestrial airglow, scattered artificial 
light, zodiacal light, scattering and emission by interstel- 
lar dust, emission by diffuse gas, and unresolved back- 
ground objects. In this section we describe a novel 
technique to measure the spectrum of scattering by the 
Galactic ISM. We use 92,000 sk y spectra fro m the Sloan 
Digital Sky Survey (SPSS. lYork et al.ll2000D . correlating 
their intensities against independently measured 100 fim 
emission to isolate the components of the DGL associated 
with interstellar dust. In a companion paper (Brandt & 
Draine 2011, in preparation) w e use the full-sky Ha map 
compiled by iFinkbeinerl t|2003l ) to isolate the component 
of the DGL associated with emission by diffuse H II. 

2.1. The SDSS Sky Fibers 

The Seventh Data Release of the SDSS 

(|Abazaiian "etani2009h contains more than 1.6 million 
spectra of stars, galaxies, and quasars, making it by far 
the largest such dataset ever assembled. The spectra 
were taken by plugging 640 fibers into a ^7 deg 2 plate, 
with each 2.96" fiber feeding the light of its target object 
into a pair of spectrographs. Each group of 640 spectra 
was then calibrated a nd sky-subtracted by the SDSS 
spectroscopic pipeline (jStoughton et all 120021 : Buries & 
Schlegel unpub.). 

To obtain an accurate measurement of the sky back- 
ground, a minimum of 32 fibers on each plate were placed 
on blank sky regions. These positions were relatively uni- 
formly distributed over each plate, and required no detec- 
tion by the photometric pipeline to 5tr in any band (Lup- 
ton, priv. commun.). A few (~1%) of the sky fibers were 
erroneously placed over bright sources, but these were 
flagged and removed by the reduction pipeline, leaving a 
total of about 92,000 blank sky spectra used to compute 
the background flux. The sky fibers were used to con- 



2 



struct a "supersky" spectrum for each plate scaled to unit 
airmass. This spectrum was then rescaled to the airmass 
at each fiber (including the sky fibers themselves) and 
subtracted from that fiber's spectrum. Each plate also 
included 8 F dwarfs as spectrophotometric standards, 8 
F sub dwarfs as reddening s tandards, and 2 hot subd- 
warfs (|Stoughton et al|[2002fl . The residual sky spectra, 
along with all of the other spectra, were flux-calibrated 
to these standards. 

The sky spectra on a given plate show modest fibcr- 
to-fiber variation, but are dominated by noise associ- 
ated with terrestrial airglow. Hidden within this noise 
are real variations due to extraterrestrial sources: zodi- 
acal light, scattering and emission by diffuse interstellar 
dust, emission by diffuse gas, and faint, unresolved back- 
ground sources. We can isolate the components asso- 
ciated with Galactic dust using the independently mea- 
sured InfraRed Astronomy Satellite (IRAS) 100_um map, 
reduced and calibrated by iSchlegel et all (|1998| ) , here- 
after SFD. The bulk of this emission comes from ther- 
mally radiating dust grains heated to ~18 K by starlight. 
By tracing interstellar dust illuminated by Milky Way 
stars, the 100 /im map allows us to measure the spec- 
trum of the Galaxy in scattered light. 

2.2. Correlating Against 100 [im Intensity 

In 1983, IRAS mapped the entire sky a t 12, 25, 60, 
and 100 microns (|Neueebauer et all 11984 ). SFD later 
smoothed the 100 /im map and corrected it for point 
sources and zodiacal light, creating a map of diffuse 
Galactic infrared emission. They f urther used COsmic 
Background Explorer (COBE) data (jBoggess et al.lll992D 
at 240 /im to estimate the temperature of the dust, ulti- 
mately producing a ~6' resolution map of 100 /im emis- 
sion and a lower resolution temperature map suitable 
to convert 100 /im emission into a dust column density. 
This map was intended mainly to estimate and correct 
for Galactic extinction. Here, we use it to correlate il- 
luminated dust with residual intensity in the SDSS sky 
fibers. 

In the optically thin limit, we expect the intensity 
of scattered starlight to be proportional to the col- 
umn density of dust times the intensity of the illumi- 
nating starlight. For dust with opacity oc v$ near 
100 /im, the intensity of the illuminating (and hence, 

scattered) starlight should be p roportional to T d ^. 
iPlanck Co llaborati on et all ()201lD have found that the 
A > 100 /im emission is well- approximated by dust with 
f3 w 1.8. We therefore expect the SDSS sky fiber resid- 
ual intensity to be roughly proportional to Tioo^m-^fust) 
where the optical depth at 100 /im, Tioo^m, is propor- 
tional to the dust column density. We expect the 100 
/im intensity itself to have a temperature dependence. 
For 18 K dust radiating at 100 /im, hv/ksT a 8 > 1, so 
that 
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(1) 



We use the measured 100 /im intensity to trace the prod- 
uct of the intensity of the illuminating starlight and the 
dust column. In practice, temperature variations are suf- 
ficiently small that using the SFD column density times 
T 5 8 produces results indistinguishable from those simply 



using 100 /im intensity. 

Two complications prevent us from assuming a linear 
relationship between sky fiber residuals and 100 /tm in- 
tensity: 

1. Such a model neglects self- absorption by optically 
thick dust. 

2. The spectroscopic pipeline has already subtracted a 
scaled sky spectrum from each fiber, which includes 
a component of the DG1. This component differs 
from plate to plate. 

We avoid the first problem by excluding spectra and en- 
tire plates where the dust is optically thick to visible 
light, with Ay > 0.5 according to SFD; our results are 
insensitive to the precise value of this threshold. For a 
dust temperature of 18 K, this corresponds to 100 /im 
emission exceeding ~10 MJysr -1 . We address the sec- 
ond point by assuming residual optical intensity to be 
proportional to the excess 100 /im intensity relative to 
the the average over that fiber's plate. Our model is then 



A/a, sky, j,p — Ct\ [(vlvJlOOf 



100 fj.m)p\ 



(2) 



where I\. s ky,j,p is the residual intensity in sky fiber j on 
plate p at wavelength A, (/„)ioo ^m,j,p is the 100 /im in- 
tensity at fiber j's location, {) p denotes an average over 
the sky fibers on plate p, and a\ is a dimensionless num- 
ber that describes the relative strength of scattered and 
thermal emission. We then solve for the best-fit spectrum 
of coefficients a\, our correlation spectrum. Defining 
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the maximum likelihood estimate for a\ is 



a 



3,P 



\V,j,P 



j,p \V,3,P 



and its variance is 



E 



3,P 



(3) 
(4) 



(5) 



(6) 



j,p Kv,3,p 



In Equations (0 and ([6]), a\ ■ is the variance of y\,j, p 
as estimated by the SDSS pipeline. 

By using the residual, sky-subtracted intensity, our 
model adopts the flux calibrations performed by the 
SDSS spectroscopic pipeline. Beginning with the Sixth 
Data Release, SDSS spectra have been flux-calibrated 
to PSF magnitudes appropriate for point sources rather 
than fiber magnitudes appropriate for extended sources. 
Dividing by the SDSS fiber a perture therefore gives an 
incorr ect intensity. Figure 4 of Adel man-McCarth v et"all 
(2008) shows the difference between PSF and fiber mag- 
nitudes of point sources to be very nearly Gaussian with 
a mean of 0.35 magnitudes. The typical seeing for spec- 
troscopy was poor (~2"); as a result, the aperture cor- 
rection was found to be a very weak function of wave- 
length. We therefore recalibrate all of our spectra to 
fiber magnitudes using the average flux conversion fac- 
tor of 10 - 35 / 2 - 5 = 1.38 before dividing by the 2.96" fiber 
aperture. 
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We would like to apply the model in Equations ©- 
(|6|) to each wavelength observed by SDSS. However, the 
wavelengths that were observed differed slightly from 
night to night. This variation is only a few tenths of 
a percent, but is sufficient to blur spectral features if not 
removed. We therefore define a new wavelength array 
of 4000 elements, similar to but slightly larger than the 
3852 element array in most SDSS spectra. We use cu- 
bic splines to interpolate all spectra and errors onto this 
new array. The interpolation introduces a correlation 
between neighboring wavelength elements. In Section 
14. 1[ we discuss this in detail, and demonstrate that our 
measured spectra and errors are statistically very well- 
behaved. 

Figure Q] illustrates our model for the spectrum of scat- 
tered light. To increase the signal-to-noise ratio in Fig- 
ure [TJ we have binned each spectrum's ^60 wavelength 
elements from 6900-7000 A. Each 2.96" patch of blank 
sky thus contributes a single point: its residual inten- 
sity averaged over the interval from 6900 to 7000 A. Af- 
ter discarding bad sky fibers flagged by the spectroscopic 
pipeline along with fibers and plates with (I^ioo/im > 10 
MJysr -1 (corresponding to Ay > 0.5), we are left with 
nearly 90,000 points representing over 5 million intensity 
measurements. We indicate these points by logarithmi- 
cally spaced contours where their density is extremely 
high. We then fit them with Equation ([2]), correlating 
the residual optical intensities against their fibers' excess 
100 /im emission. There is a positive correlation between 
the two quantities (a\ > 0), driven by the entire ensem- 
ble of points and significant at more than 70 a. 

The best-fit slope in Figure [T] is biased low by its ne- 
glect of errors in the 100 /zm map and by structure un- 
resolved by IRAS, but by a factor independent of wave- 
length. Chromatic effe cts due to seeing and atmospheri c 
dispersion are small ([Adclman-Mc Carthv et all [2008). 
We discuss the calibration of the correlation spectrum 
in Sections B~2l and IPl 

2.3. Sky Coverage 

The SDSS focused much of its spectroscopy on the 
Northern Galactic Cap because of that region's low lev- 
els of extinction. For the same reason, many of its plates 
are not very useful for inferring the spectrum of scattered 
light. The part of the SDSS footprint that is most useful 
for our purposes has significant small-scale variation in 
100 (im emission, and is disproportionately in regions of 
low Galactic latitude. 

Figure [2] shows the sky coverage of the SDSS spec- 
troscopic survey (upper panel), and the sky coverage 
weighted by its power to affect the fit described by Equa- 
tion ([2]) (lower panel) . For any plate, the effect on the fit 
will be roughly proportional to the variance of the 100 
/zm intensity over that plate, 

To obtain the weighted map, we have multiplied the fiber 
density by the ratio of crl 00fim p to its average over the 

entire sample, (ofou/mi)- 

Stripe 82 is clearly visible in the lower-left (and ex- 
treme lower-right) parts of each map. This region, de- 
fined by -50° < RA < 59°, -1.25° < 5 < 1.25°, 
was repeatedly imaged as part of the SDSS Supernova 



- 6900-7000 A s ^pL** 
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Fig. 1. — Scatter plot showing the correlation between 100 (im 
intensity relative to the mean over an SDSS sky fiber's plate and the 
residual sky intensity averaged over the interval from 6900 to 7000 
A (Equation 0. We use logarithmically spaced contours where 
the density of points is high. The data demand a non-zero slope, a 
component of the sky background associated with interstellar dust, 
at more than 70 a. 

Surve y (lAdelman-McCarthy et al.l 120061 : iFrieman et al.l 
120081: lAbazaiian et al.l 120091) . It was also the target of 
sever al "special" plates not part of th e main SDSS sur- 
veys (jAdelman-McCarthv et al.l 120061 ). Stripe 82, along 
with a variety of patches at modest Galactic latitude 
and a swath at the southern limit of SDSS's main foot- 
print (in Galactic coordinates), contributes the bulk of 
our measured signal. 

3. RESULTS: SPECTRA 

Here we present our correlation spectra of the DGL, 
a\, computed as described in Section[2j We present three 
sets of spectra: the continuum of the scattered spectrum, 
and the spectra from 4830-5040 A and from 6530-6770 A 
to show emission lines. We also divide the sky into three 
regions by Galactic latitude and longitude to explore the 
spatial variation of the scattered light. 

3.1. Continuum of the Correlation Spectrum 

We obtain the spectrum of the DGL associated with 
100 /im emission using the method described in Section 
[272] By fitting Equation $2$) to the SDSS residual sky 
spectra, we obtain a dimcnsionlcss coefficient a\ at each 
wavelength relating optical intensity to 100 /im emission. 
We then mask the nebular emission lines Ha, H/3, [N n] 
A6550, [Nil] A6585, [Sn] A6718, and [Sn] A6733 and bin 
the correlation spectra in intervals of 50 A to increase 
the signal-to-noise ratio. We do not mask auroral lines 
like [O i] A6300, which are uncorrelated with the 100 /im 
intensity. As we show in Section [4. 11 our errors are inde- 
pendent except for an interpolation effect, and they are 
normally distributed. The errors on the binned spectra 
pass the same statistical test. 
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<1 2 4 8 16 32 64 128 + 

Weighted Sky Fiber Density [deg~ 2 ] 

Fig. 2. — Top panel: the sky coverage of the SDSS blank sky 
fibers in Galactic coordinates, centered around (l,b) = (0,0). Bot- 
tom panel: the sky coverage of the blank sky fibers weighted by 
f7 l00Mm, P /( cr iOO M m) (Equation 0), which provides an approxi- 
mate measure of each fiber's influence on the correlation spectrum. 
The bulk of the SDSS footprint is near the North Galactic Pole, 
with additional concentrated sampling in Stripe 82 near the South 
Galactic Pole. The useful survey footprint is much more heavily 
weighted toward equatorial regions. 

Figure [3] shows the continuum spectrum of the DGL 
computed using all the sky fibers. The spectrum is very 
blue, yet it shows a clear 4000 A break characteristic of 
old stellar populations. Broad Mg and Fe b absorption 
are also visible just blueward of 5200 A. All of these char- 
acteristics are consistent with a continuum of scattered 
starlight. A simplified radiative transfer calculation, dis- 
cussed in Section 15.11 confirms scattering as the source 
of the DGL and shows that it may be used to discrimi- 
nate between dust models. The four plotted curves use 
two estimates of the continuum of the inter stellar radia- 
tion fi eld (I SRF) and two dust models, t he IZubko et all 
(|2004f) and iWeingartner fe Draind (|200lD models (here- 
after ZDA04 and WD01); an excess of large grains in the 
latter produces a redder scattered spectrum. As we dis- 
cuss in Section r4.2l errors and small-scale structure in the 
100 /im map bias our recovered correlation spectrum low 
by an unknown factor, which we estimate to be 2.1 ± 0.4 
(see Section |4~2]) . We have therefore scaled the radiative 
transfer models to a common level to show their different 
shapes. The scale factors are indicated on the plot, and 
are all consistent with our estimate for the bias. 

Figure 2] shows the continuum spectra from restricted 
areas of the sky. The left panel divides the sky by Galac- 
tic latitude, with the ranges chosen to have comparable 
100 /im emission summed over our sky fibers. In the high- 
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Fig. 3. — The full-sky continuum of the scattered DGL, calcu- 
lated as described in Section [2] We have masked nebular emission 
lines and binned the remaining spectrum into 50 A intervals. The 
continuum is very blue and shows a clear 4000 A break and Mg and 
Fe b absorption just blueward of 5200 A; all of these features sug- 
gest scattered starlig ht. Simplified radiative transfer calculations, 
described in Section 15.11 confirm this and demonstrate that the 
DGL can discriminate between dust models. Because our correla- 
tion spectrum is biased low by an unknown factor (Section I4.2jl . 
we have scaled all spectra to a common level. The dashed lines 
use a featureless ISRF, while the solid lines use stellar population 
synthesis models. The ZDA04 model has fewer large grains and 
gives a better fit than WD01. We find no evidence of Extended 
Red Emission (ERE) around 7000 A ( Section [5~4l . 



est latitude bin, this flux is distributed over a much larger 
number of sky fibers and the resulting measurement is 
considerably noisier (see Figure [5). Still, the errors are 
reliable (see Section l4~Tj) . and the shape differs from its 
values at lower latitude by many sigma. It is not clear 
whether the spatial variation is due to variation in the 
dust properties, in the illuminating starlight, systematic 
effects that become important when the signal is weak, 
or some combinatio n of these factors. In particular, 
lYahata et al.1 (|2007t ) find a correlation of SDSS galaxy 
density with SFD 100 /im emission in regions of low ex- 
tinction, consi stent with extragala ctic contamination of 
the SFD map. lYahata et al.l (|2007l ) report that this con- 
tamination corresponds to an inferred Ay « 0.01, or 
-Tioo/im ~ 0.2 MJysr -1 . This would represent about 15% 
of the 100 /im intensity in our sky fibers with |6| > 50°, 
and an even larger fraction of the intraplate variation to 
which we are sensitive. 

The right panel of Figure HI shows the continuum spec- 
tra for different regions in Galactic longitude. We take 
the longitude I to run from —180° to 180°, and choose 
our regions to be equal in area and symmetric about the 
Galactic center. About half of our signal comes from the 
region opposite the center, with 120° < \l\ < 180°. As 
with Galactic latitude, there are significant spatial vari- 
ations in the correlation spectrum. It is unclear whether 
the redder spectrum in the direction of the Galactic cen- 
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Fig. 4. — The continuum spectra of the scattering component of the DGL for different regions of sky. These are identical to Figure[3]for 
restricted ranges of I and b; note that we take —180° < / < 180°. The continuum shows some variation with Galactic latitu de (left panel) 
and looks qualitatively different at the highest latitudes, potentially due to extragalactic contamination of the 100 (im map 1 Yahata et al l 
120071) . The spectrum is also significantly redder toward the Galactic center (right panel). This spatial variation could be due to differences 
in the dust composition, in the illuminating starlight, or uncorrected systematic effects. 



ter is due to larger grains, redder illuminating starlight, 
or some other effect. This spectrum is the most sen- 
sitive to the cutoff at high optical depth (Section 12. 2[) , 
which may indicate significantly reddened illuminating 
starlight. 

3.2. Emission Lines 

Figure [5] shows the correlation spectrum in the wave- 
length ranges 4830-5040 A and 6530-6770 A, computed 
using all of the sky fibers without binning the a\. The 
spectrum of the DGL exhibits strong nebular emission 
lines which we masked to show the continuum in Fig- 
ures E] and H We detect Ha, H/3, [Nil] A6550, [Nil] 
A6585, [S n] A6718, and [S n] A6733 at high significance, 
but see only weak emission in the [O in] A5008 line ex- 
cited by early O-type stars. Unfortunately, the [O n] 
A3727 line lies just outside the SDSS wavelength range. 
These emission lines likely represent scattered photons 
that were originally emitted from H II regions and the 
diffuse warm ionized medium (WIM). 

We list line strengths in TableQ]measured as equivalent 
widths, which are unaffected by bias factors fSection l4.2[) . 
We use an interval of 6 SDSS wavelength elements, cor- 
responding to a velocity range of 400 kms -1 , to measure 
line intensities. For the lines between 6500 and 6800 A, 
we use the average intensity over the interval from 6600 
to 6700 A as our estimate of the continuum. For H/3, 
we use 30 wavelength elements on each side of the line, 
running from 4829-4858 A and from 4864-4893 A. The 
ratio of the [Nil] doublet, [Nil] A6585 to [Nil] A6550, 
provides a check on our recovered spectrum; our mea- 
sured value of 2.8 ± 0.6 matches the ratio of 3 expected 
from the Einstein A coefficients. 

In order to measure the strengths of the Balmer lines 



in emission, we need to account for the fact that they 
appear in absorption in stellar spectra. Fortunately, the 
strength of the Balmer absorption lines in composite stel- 
lar spectra is closely correlated with that of the 4000 A 
Calcium break. We measure the 4000 A break m our con- 
tinuu m spectrum and use a ran ge of model stellar spectra 
from lBruzual & Charlotl (|2003| ). hereafter BC03, to fit a 
linear relationship between the strength of the break and 
the Balmer equivalent widths. We use models of 6 Gyr 
of constant star formation, single stellar populations of 
2.5 Gyr, 5 Gyr, and 11 Gyr, and two exponential star 
formation histories, each with metallicities of 0.02 and 
0.008 (Zq and 0.4Z Q ); a combination of these models 
should provide a reasonable fit to stellar populations in 
the Solar neighborhood. Defining J4000 to be the ratio 
of the integrated intensity between 3850 and 4000 A to 
the integrated intensity between 4000 and 4150 A, and 
computing the continua and the line widths as described 
above, we find best-fit relationships of 

^ EW|H^) \ ^ _ 2 2§4ooq + o 17 and (8) 

(^^l)^-1.55 4 ooo-0.19 (9) 

for the model stellar spectra. The root-mean-square scat- 
ters of the BC03 equivalent widths around the fits given 
by Equations ((U) and ((9|) are ~0.06 A and -0.1 A, re- 
spectively, significantly smaller than the errors in our 
measurements of the DGL. Correcting for stellar absorp- 
tion increases our measured Ha and H/3 line strengths 
by -20%. 

The corrected equivalent width of Ha in emission is 
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Fig. 5.— Emission lines in the DGL. H/3, Ha, [Nil] A6550, [N II] A6585, [S II] A6718, and [S II] A6733 are all detected with high significance, 
while we find little evidence of emission in [O III] A5008. The lack of [O III] indicates a relatively soft spectrum of ionizing photons compared 
to average Hn regions, while strong [S II] and [Nil] indicate warm, ~8,000 K gas. 



TABLE 1 



Equivalent Widths and Line Ratios in the DGL 




Line 


Equivalent Width [A] 


Energy Ratio of Line to 






Ho 


H/3 


H/3 A4863 


4.8±0.7 a 


0.38 ±0.05 


1 




Om] A4960 


1.6 ±0.6 


0.14 ±0.05 


0.38 ±0.15 




Om] A5008 


0.8 ±0.6 


0.07 ±0.05 


0.19 ±0.15 




Hci] A5877 


0.3 ±0.8 


0.03 ±0.07 


0.09 ±0.21 




N ill A6550 


2.4 ±0.5 


0.19 ±0.04 


0.50 ±0.12 


Ha A6565 


12.5±0.5 a 


1 


2.64 ±0.38 


[No] A6585 


6.6 ±0.5 


0.53 ±0.04 


1.40 ±0.22 


[S n] A6718 


5.7 ±0.4 


0.45 ± 0.04 


1.20 ±0.19 


[S ii] A6733 


4.3 ±0.4 


0.34 ±0.04 


0.91 ±0.15 



TABLE 2 
The Eight Nearest O Stars' 1 



a Corrected for stellar absorption using Equations (01 and j9) . 

consistent wi th its value of 11 A in the local ISRF 
(jDraind I2011L Table 12.1), obtained by integr ating; the 
full-sky Ha map compiled by Finkb einer (2003) fro m th e 
Wisconsin Ha Mapper (WHAM. iRevnolds et all 12 002). 
South ern Ha Sky Survey Atlas (SHASSA. lGaustad et alj 



20011), and Virginia Tech Spectral-Line Survey (VTSS, 
Dennison et al.lTi998l) . The line ratios provide a probe 



of the average physical conditions in the local ISM. The 
strength of the singly ionized [N n] and [S n] lines and 
weakness of [Om]A5008 and [Hel]A5877 indicate that 
most of the S, N, and O are singly ionized, while the He 
is largely neutral. This implies a lack of photons with 
hv > 24.6 eV in the Solar neighborhood, which can be 
understood from the lack of stars of spectral type 08 
and earlier within 300 pc of the Sun. The nearest eight 
O stars are listed in Table [2] We discuss the physical 
conditions of the local ISM in more detail in Section [SJ 
Figure [B] shows the strength of the emission lines Ha, 
[Nil] A6550, [Nil] A6585, [S n] A6718, and [S n] A6733 
relative to 100 /im emission for different ranges of Galac- 



GOS ID 



Other ID Spectral Type Distance (pc) 



G006.28±23.5901 
G203.86-17.7401 
G206.45-16.5901 
G202.94+02.2001 
G255. 98-04.7101 
G262.80-07.6901 
G195. 05-12.0001 A Ori A 
G206. 82-17.3401 a Ori AB 



C Oph 
5 Ori A 
C Ori A 
15 Mon 

C Pup 
7 2 Vel 



09.5V 
09.5V 
09.71b 
07V 
041 
09:1: 
08III 
09.5V 



221±25 
239tg 
309±f° 

349±|| 

o 61 +89 

380 tlf 



a References: Mai'z-Apcllaniz ct al. (2004); Mai'z Apcllaniz et al. 
(2008) 



tic longitude and latitude. The lines are somewhat 
stronger at low Galactic latitude and in the direction of 
the Galactic center, though this could reflect the relative 
number of nearby H II regions rather than the physical 
conditions in the ISM. An important caveat is that, be- 
cause of spatial variations in the 100 /^m intensity, the 
correlation spectra in different regions of the sky need not 
share the same bias factor (Sections 14.21 and 14.31) . The 
ratios of the [S n] and [N n] lines to Ha, which are robust 
to calibration difficulties, vary little across the sky. 

4. CALIBRATION AND MEASUREMENT ERRORS 

The standard errors on our correlation spectra oc\ are 
very nearly normally distributed with the variance given 
by our maximum likelihood estimator, up to a constant 
factor due to our interpolation of the original SDSS spec- 
tra onto a common wavelength array. The absolute cal- 
ibration of the a\ is more problematic. Our neglect 
of (unknown) measurement errors in the 100 /jm inten- 
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Fig. 6. — Emission line spectra of the dust component of the DGL as functions of Galactic longitude and latitude. The line ratios and 
equivalent widths depend only weakly on position, though the equivalent widths of all emission lines seem to decrease away from the 
Galactic center (green line, right panel). 



sity and structure unresolved by IRAS introduces a bias, 
nearly independent of wavelength, which we estimate to 
be a factor of 2.1 ± 0.4. We demonstrate both of these 
results below. 

4.1. Measurement Errors 

The errors in our correlation spectra are derived from 
fits to about 90,000 intensities over the full sky. The sky 
spectra are first interpolated onto a common wavelength 
array, introducing a correlation between the intensities 
at neighboring wavelengths. The difference between the 
intensities I\ at wavelength elements i and i + 1 ranges 
from i to 1 time(s) its value without interpolation (a 
factor of i corresponds to the original and final arrays 
being offset by half a wavelength increment). Because the 
offsets between the original and resampled wavelengths 
are random, the denominator of this fraction is uniformly 
distributed between 1 and 2, which reduces the average 
difference (li+i — h) by a factor of 



'i^orig 



{I 



i+1 



H Jintcrp 



— = In 2 

t 



(10) 



This correlation will propagate through the individual in- 
terpolated spectra to the ct\. If the element-to-element 
variations in ct\ are dominated by Gaussian measure- 
ment errors, we expect the quantities 



H+l 



OU 



(In 2) 



(11) 



i+l 



to be normally distributed with unit variance; on and 
are given by Equations ([3]) and ©, respectively. 

Figure [7] shows that the distribution of Ao^ as defined 
by Equation is exceedingly well-fit by a normal dis- 
tribution with zero mean and unit variance. We have 



masked 67 of the 3999 normalized intensity differences 
which lie near the emission lines Ha, H/3, [Nil] A6550, 
[N n] A6585, [S n] A6718, and [S n] A6733. At these wave- 
lengths, real spectral features contribute to the element- 
to-element variation in ot\. Even with 3932 values of Aq^ 
across the wavelength range, a Kolmogorov-Smirnov (K- 
S) test only detects deviation from Gaussianity with 83% 
confidence. 

The exceptional agreement shown in Figure [7] gives us 
confidence that our measured errors are reliable and in- 
dependent (up to a factor of In 2 from interpolating). We 
have verified that coadding neighboring wavelength ele- 
ments, as we did to smooth the continuum in Section [3~T 
and to compute the equivalent widths of lines in Section 
does not affect the statistical properties of the errors. 



4.2. Biases in the Correlation Spectra 

Our model of the correlation between 100 /im intensity 
and optical intensity (Equation (|2"|)) includes errors in the 
SDSS sky fiber residuals, but neglects (unknown) errors 
in the 100 /im intensity at each fiber's location and struc- 
ture unresolved by IRAS. This introduces a wavelength- 
independent bias to our recovered spectra (Figures [3] - 
[5]). We first derive an expression for the bias and then 
estimate its value. 

The correlation spectrum ot\ (Equation ([2"])) is derived 
from a x 2 minimization. The maximum likelihood values 
of the a\ and their variances are given by Equations 
and (JB]), with Xj tP and yx,j,p defined in Equations © and 
((31) ; Xj.p is the excess 100 /im intensity in fiber j relative 
to the average on plate p, yx,j,p i s the residual sky fiber 
intensity XI\ at wavelength A in fiber j on plate p, and 
°a j p i s ^ s variance as estimated by the SDSS pipeline. 
We let £j tP denote the true excess 100 /mi emission at 
sky fiber j on plate p, so that the measurement error, 
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Fig. 7. — Normalized diffe ren ce between «a at successive wave- 
length elements (Equation mil l) in the full sky correlation spec- 
trum, masking wavelengths near emission lines. The factor of In 2 
arises from the fact that we interpolate all spectra onto a common 
wavelength range, which introduces a co rrela tion between neigh- 
boring wavelength elements (see Section 14.11 1. There is no other 
scaling. The resulting distribution of 3932 points is very nearly 
Gaussian; a K-S test rejects the consistency of the two distribu- 
tions with just 83% confidence. 

including the effects of unresolved structure, is 8j tP = 

x j,p ~ £j,p- _ _ 

Assuming the model in Equation (J2J) to be correct, we 

may write yx,j,p as ax£,j. P + £\,j,p, with Sx,j, P represent- 
ing the measurement error in the sky fiber intensity. We 
further assume the error terms £x,j,p and 6j iP to be mi- 
correlated with zero mean and invoke the large number 
of sky fibers (nearly 10 5 ) to neglect sums of 6 and e. The 
first factor in Equation ([5]) becomes 



E 

3,P 



{yx,j,p) ( x j,p) 



a 



ax 



hp 



X,y,j,p 



E CT 2 

j,p *,y,3,p 



(12) 



This is the same value we would measure with b^ p = 0. 
Noting that the mean 100 (im excess, (xj_ p ), is zero by 
construction, our estimate of ax is biased low by a factor 

-l 



ax(x,y) 
ax(£,,y) 




(13) 

-(J 2 ) 



where cr 2 is the observed variance in x and cr| 
would be its value with no measurement errors. Thus, 

ax(x,y) 



l 



rr2 



(14) 



Note that this bias is independent of wavelength and of 
the errors in the sky fiber intensities. We have empiri- 
cally verified the latter by adding noise to the intensities; 
we recover the same correlation spectra to within the er- 
rors. 



4.3. Calibrating the Correlation Spectra 

An unbiased estimator for ax, for example a likelihood 
function of the form 



C(xj,yxj\ax) = J C(xj\x) C(yx,j\x,ax) C[x) dx , (15) 

would remove all calibration issues. The SDSS pipeline 
provides an estimate of C(yx,j\x,ax), the likelihood of 
measuring sky fiber residual yx in sky fiber j given a 
correlation spectrum ax and true 100 /jm intensity x. 
Unfortunately, we have almost no information on the er- 
rors in the 100 /an intensity to estimate C(xj\x), the 
likelihood of measuring excess 100 /im intensity Xj given 
a true value x, and can only guess at a prior, C(x), from 
the 100 /xm map itself. 

Because of these difficulties, we use two alternative and 
independent approaches. We first construct an estimator 
that we expect to be asymptotically unbiased, assuming 
the error in residual 100 /im intensity to depend weakly 
on the true residual value at a fiber position. This is a 
reasonable assumption, particularly because we are sub- 
tracting the mean 100 /im emission over a plate; an re- 
value of zero does not correspond to zero intensity. We 
then show the results of our radiative transfer calcula- 
tions assuming ZDA04 dust and a plane-parallel galaxy. 

Under the assumption of uniform errors in the residual 
100 /mi intensity, Equation (fT4]) shows that restricting 
the sample to fibers with large |a;| (and therefore large 
<7 2 ) will give an asymptotically unbiased estimate. We 
therefore recalculate a\ between 6600 and 6700 A using 
only the fibers with \x\ > |ac| m m. Figure [8] shows a\ as a 
function of |x| m i n , calculated by varying |x| m i n from to 
4 MJysr ; as expected, a\ increases with |x| m j n - Figure 
[8] suggests a bias factor of at least ~1.7— 2. This is also 
supported by our radiative transfer model, discussed in 
detail in Section ISTTl which suggests a bias of wl.7— 2.4. 
We conservatively adopt a bias factor of 2.1 ± 0.4, indi- 
cated by the shaded region of Figure [51 to calibrate our 
correlation spectrum. 

5. DISCUSSION 
5.1. Scattered Light and Dust Models 

Many features of the DGL - the 4000 A break, broad 
Mg + Fe b absorption, and a much bluer continuum 
than that of stars with these spectral features - support 
the hypothesis that the DGL is dominated by scattered 
starlight. This is hardly surprising, as the spectrum was 
derived by correlating residual optical intensity with 100 
/im intensity over small spatial scales. A simplified ra- 
diative transfer calculation confirms that our spectrum is 
consistent with dust scattering and allows us to use the 
DGL to discriminate between dust models. 

Our model of the DGL uses an infinite plane-parallel 
galaxy with a Gaussian vertical distribution of dust, 



p dust cx exp 



z 

~2^ 



(16) 



with a = 250 pc (|Malhotral 119951: iNakanishi fe Sofud 
2003). We approximate the stellar distribution as 
the sum of two exponentia l distributions with scale 
heights of 300 pc and 1350 pc (|Binnev fe MerrifieldlU998l: 
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Minimum |I 100(im - <I 100fim >| [MJySr" 1 ] 

Fig. 8. — The correlation spectrum ct\ between 6600 and 6700 
A using only points with large excess 100 fim intensities. The red 
curve indicates the number of sky fibers satisfying each 100 fim 
cutoff, ff the errors on residual 100 /im intensity depend weakly on 
fiber position, this estimator will be asymptotically unbiased. We 
overplot several single-scatterin g ra diative transfer models, which 
we describe in detail in Section 15.11 All assume ZDA04 dust, but 
with different dust densities and different fSRFs. We conservatively 
adopt the bias factor of 2.1 ± 0.4 indicated by the orange shading 
to calibrate our correlation spectrum. 

IGilmore fc Reidl I1983H . The 300 pc component domi- 
nates the distribution with about 90% of the stars. 
We use two estimates of the stellar emission spectrum: 

1. A model th at reproduces the local ISRF of 

iMathis et all (|1983h . hereafter MMP83, in the mid- 
plane, and 

2. A stellar population synthesis model from BC03, 
with solar metallicity and an exponential star for- 
mation history over 12 Gyr. 

The MMP83 ISRF for A > 2460 A is approximated as 
a sum of three dilute blackbodies of T = 3000, 4000, 
and 7000 K, with dilution coefficients W = 5 x 10~ 13 , 
1.65 x 10 -13 , and 10 -14 , respectively. We comput e th e 
attenuation of the stellar emission using Equation (IA3|) . 
setting the attenuated spectrum equal to the local ISRF. 
In this way, we "de-redden" the MMP83 ISRF to obtain 
the stellar source spectrum. 

We show all of the ISRF spectra in Figure [SI The peak 
of the de-reddened MMP83 \F X is at about 8000 A, cor- 
responding to a temperature of ^4600 K. The spectrum 
is similar to a BC03 model with a 5 Gyr star forma- 
tion timescale, though MMP83 has less UV emission. 
This is probably a result of the relative lack of early 
O stars in the Solar neighborhood and of our model's 
neglect of extinction from d ust in a young star's birth 
cloud (jCharlot fc Falll l2000h . The UV discrepancy be- 
comes much more serious for BC03 models with more 
extended star formation histories. This UV excess sig- 
nificantly increases the dust heating and decreases the 
ratio of scattering in the optical to emission in the far- 
infrared. For a star formation timescale of 5 Gyr, this 
is a ~30-40% effect relative to the de-reddened MMP83 
model (see Figures 151 and [TU)) . 

Theoretical and observational estimates of the lo- 
cal star formation history favor roughly constant star 
formation rates over ~10 Gyr (He rnandez et all 120011 : 



MMP83 ISRF 

De-Reddened MMP83, t vo = 0.05 

De-Reddened MMP83, t V0 = 0.15 

BC03, Z , SFR « exp(-t/9 Gyr) 

BC03, Z , SFR « exp(-t/5 Gyr) 




1000 10000 
Wavelength [A] 

Fig. 9. — Comparison of the MMP83 ISRF, de-reddened using 
our radiative transfer models, and BC03 models with solar metal- 
licity and exponential star formation histories over 12 Gyr. The 
BC03 models have significant excesses in the UV that increase dust 
heating and reduce the ratio of the optical scattering to far-infrared 
emission. 

ICignoni et al.l r2006). Such models do not agree with our 
measured spectrum of the DGL and, using our radiative 
transfer model, would produce a very different ISRF from 
MMP83. This may indicate that a substantial fraction 
of the illuminating starlight originates relatively far from 
the Solar neighborhood, it may be a result of our simpli- 
fied radiative transfer, or it may indicate that stars in the 
Solar neighborhood are generally older than is currently 
thought. The spatial variation of the DGL (Figure [4j 
shows that the geometry of the radiative transfer prob- 
lem is likely to be important. A more detailed model 
of the Galaxy could better constrain the average stellar 
source spectrum. 

Our radiative transfer calculations neglect multiple 
scatterings. This is a good approximation at high Galac- 
tic latitude where the optical depths are low. It becomes 
poor near the midplane, but SDSS has very little sky 
coverage near the Galactic plane (Figure [2]) . We assume 
all absorbed starlight to be reradiated isotropically in 
the infrared and use a Henyey-Greenstein phase func- 
ti on for scatter i ng in the optical. We use the dust model 
of iDraine &: Lil ([20071 ) to convert total infrared power to 
IRAS 100 /im bandpass power, with 

(^Wm = (0-52±0.05)/ T iR . (17) 

The central value corresponds to their model with an in- 
cident starlight intensity 80% of th at in the Solar neigh- 
borhood, which IDraine fc Ll ([2007D found to provide the 
be st fit to the average far -infrared spectrum measured 
by iFinkbeiner et al.l JT999). The confidence interval in 
Equation (|17l) includes models from 0.5 to 1.5 times the 
local starlight intensity. Once a dust model is specified, 
with wavelength-dependent cross-sections, albedos, and 
anisotropy parameters, our model for the scattered light 
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Fig. 10. — Single-scattering radiative transfer calculations assum- 
ing a plane-parallel galaxy. The results are much more sensitive 
to the assumed dust model than to the parameters of the galaxy 
model. The WD01 model produces more scattering at longer wave- 
lengths than the ZDA04 model because of an excess of large dust 
grains. 



spectrum has no free parameters. We derive the relevant 
equations in the Appendix and evaluate them numeri- 
cally. 

Figure [10] shows the results of our calculations to be 
remarkably insensitive to the details of the galaxy mod- 
eling (other than the assumed stellar source spectrum). 
Though not shown, models with a single exponential 
distribution of stars and with an exponential, rather 
than Gaussian, dust distribution, are nearly indistin- 
guishable from the present models. For esc \b\ > 1.4 
(\b\ < 45°) the H i has iV H « 2.9 x 10 20 csc|fe| cm~ 2 
(|Dickev et al.l [l978() falling below this relation for b > 
45°. For E(B - V)/N R = 5.8 x 10 21 cm 2 mag- 1 and 
Ry = Ay/E{B — V) = 3.1, we would then expect 
Ty = 0.17 esc |6| for \b\ < 45°, and lower values at 
b > 45°. We use Ty = 0.15 esc |6| for our fiducial dis- 
tribution. 

Figure ITT1 shows that the optical depths at the SDSS 
sky fibers, computed using the SFD estimates of E(B — 
V) and assuming Ry — 3.1 dust, are roughly bracketed 
by ry = 0.05 esc |6| and Ty = 0.15 esc |6|. This suggests 
that the sky fiber locations were slightly biased toward 
lower-than-average H I column densities. 

While the details of the model galaxy have little effect 
on the predicted spectrum, the dust model matters a 
great deal. The WD01 model appears to have too many 
large grains, giving too much scattering at long wave- 
lengths. The size distribution of ZDA04 brings our model 
into much better agreement with the data (see Figure [3]). 
While our model predicts the scattering spectrum to be a 
very weak function of Galactic latitude, variations in the 
dust properties, geometry, or illuminating starlight could 
result in much larger differences. Possible spatial varia- 
tions in 100 /jm errors and small-scale structure (Section 



Fig. 11. — Optical depths at each sky fiber computed using the 
SFD estimates of E(B — V) and assuming Ry = 3.1 dust. The 
dust distribution is roughly bracketed by ry = 0.05 esc |b| and Ty = 
0.15 esc |6|; these two models predict nearly the same spectrum of 
scattered light. 



14. 2p further complicate the interpretation of the spatial 
variations seen in Figure |4l 

5.2. Line Emission 

On the evidence presented above, we can be confident 
that the continuum of the correlation spectrum a\ con- 
sists primarily of scattered starlight. It is more difficult 
to show that the line emission is scattered rather than 
from ionized gas physically associated with the dust. The 
strongest piece of evidence is that the equivalent width of 
Ha in the correlation spectrum, 12.5 ±0.5 A, matches its 
measured value of 11 A in the local ISRF (Section 13. 2p . 
Because scattering will preserve the equivalent width of 
Ha in the ISRF incident on the dust, a stronger Ha line 
would have indicated an additional, nearly continuum- 
free component seen in direct emission and correlating 
with 100 /Ltm intensity. 

If the emission lines observed in Section 13.21 are ob- 
served mostly or entirely in reflection, they provide a 
probe of the average physics of the nearby ISM. As dis- 
cussed in Section 13 21 the strength of the [N n] and [S n] 
lines, combined with the weakness of the [O in] lines, in- 
dicate relatively few photons with hv > 24.6 eV. This 
is probably due to the lack of early-type O stars in 
the Solar neighborhood. The observed strengths of the 
collisionally-excited lines relative to recombination lines 
of similar wavelength are: 



[S n] (6718 + 6733) 

Ha 6565 
[Nn](6550 + 6585) 

Ha 6565 
[O in] (4960 + 5008) 
H/3 4863 



0.80 + 0.05, (18) 
0.72 + 0.06 , and (19) 
0.58 + 0.21. (20) 
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Fig. 12.— Upper panel: [Nn](6550+6585)/Ho vs. T. Curves 
are predicted values for indicated values of Nil/Ho relative to 
(N/H) . If 0.7 <(Nn/Hn)/(O/N) < 1, then the observed 
line ratio 0.72 ± 0.06 requires 7180 < T < 8400K. Middle panel: 
[S n] (6718+6733) /Ha vs. electron temperature T. Curves are pre- 
dicted values for indicated values of Sii/Hii relative to (S/H)q. 
The observed line ratio 0.80±0.05 and 0.7 <(S n/H Il)/(O/H) < 1 
limits the temperature to 7180 < T < 8330K. Lower panel: 
Om(4960+5008)/H / 3 vs. T. Curves are predicted values for 
indicated values of Om/Hll relative to (0/H)q. The ob- 
served line ratio 0.58 ± 0.21 and the allowed temperature range 
7180 < T < 8330K implies that Oin/Hn is between 0.039 x and 
0.15x(O/H)q. The dot in each panel indicates our best esti- 
mate: T = 7700K, Nn/H 0=0.85 (N/H) , S o/H n=O.83(S/H) , 
and Oiii/Hii=O.O8(O/H) . 



These allow us to estimate the temperature of the ISM 
where the lines originate and the state of ionization 
of S, N, and O. Figure [12] shows predicted line ratios 
as functions of electron temperat ure T, calcula ted us- 
ing Ha and H/3 emissivities from iDraind (1201 If), colli- 
sion stre ngths for N II from Hudson fc B ell (|2005l) . for 
Sll fromlTaval fc Zatsarinnvl (|201Q[). and for O in from 



10^ cnT 



lAggarwal k, Keerianl ([1999). A density 
was assumed; the results are insensitive to n e provided 
n e < 10 3 cm -3 . 

Nitrogen is not depleted in the ISM, and the first and 
second ionization potentials (14.0 and 29.6 eV) lead us 
to expect Nil/Hn to be close to the solar abundance 
(N/H)© = 7.4 x 10" 5 (lAsplund et alJ[20?rl . If we as- 
sume 0.7 <(Nii/Hii)/(N/H) < 1 we see from Figure 
[T3]that the observed [N n] (6550+6565) /Ha ratio allows 
only temperatures 7180 < T < 8400 K. 

Sulfur is not expected to be depleted in Hll regions 
or the diffuse ISM. With an ionization potential of 23.38 



eV for S ii — >S III, S n will be the dominant ionization 
stage in H II regions where He is neutral. If we assume 
S n/H II to be between 0.7 and l.Ox the solar abundance 
(S/H) = 1.45 x 10" 5 (jAsplund et al.ll2009l) . then we see 
from Figure [H that the observed [S li(6718+6733)]/Ha 
requires 7180 < T < 8330 K. 

Oxygen is only slightly depleted in the diffuse ISM, 
with ~20% of the O resident in silicates. The second 
ionization potential of oxygen is 35.1 eV, and therefore 
O in will be present only when He is ionized. A star of 
spectral type 08 or earlier is required for the He ioniza- 
tion zone to account for m ore than 50% of the mass in the 
H II region (Draine 2011). As seen in Table[2j the nearest 
such stars are 15 Mon (07V, d = 309±|| pc) and C Pup 

(041, d = 335±u P c )- The reflected light in the DGL is 
expected to originate mainly within a few hundred pc of 
the Sun, and therefore the contribution from H II regions 
should be dominated by H II regions where He is neutral 
and O is singly ionized. Figure [T3] shows that the ob- 
served strength of [O m](4960+5008)/H/3 is consistent 
with Oin/ Hn between 039 a nd 0.15 of (O/H) = 
5.4 x 10" 4 (|Asplund et all [20091 ) . The observed [Om] 
emission can be reproduced by emission from H II regions 
with T w 7700 K and O m/H II w 0.08 x (O/H) . 

He I 5877 is not detected, with a 3-cr upper limit 
He I 5877/Ha < 0.24; for T w 7700 K, this corresponds 
to Heii/Hn < 0.48, which is not a useful constraint. 

The emission lines measured in the correlation spec- 
trum a\ are therefore consistent with the line ratios ex- 
pected for H II regions within ^400 pc of the Sun, for 
electron temperature T w 7700 ± 250 K. These tem- 
peratures and ionization state s are also consistent with 
the WIM (jMadsen et al.ll2006| ). which we expect to con- 
tribute a significant fraction of the line emission, partic- 
ularly at high latitudes. 

5.3. The Fraction of Scattered Ha at High Latitude 

The previous section argued that the Ha in our corre- 
lation spectrum is scattered. We can then use the cali- 
brated correlation spectra to estimate the fraction of the 
observed Ha at high Galactic latitude that is scattered 
light. We simply integrate the correlation spectrum over 
the Ha line and multiply by the ratio of the total 100 
/mi to Ha emission over the high latitude sky fibers. 

From Figures [8] and [101 we take the ratio 
A7 A (6600 k)/vl v (100 pm) to be 0.38 ±0.07. This contin- 
uum value allows us to convert an Ha equivalent width 
in our correlation spectrum into an intensity relative to 
the 100 /im intensity. We then have 



/(Ha) 
R 



(0.090±0.017)x 



/ W A (Ha) \ / /, y (100 M m) \ 

V 12.5 A J X V MJysi- 1 J ' 

(21) 

where the Ha intensity /(Ha) is in Rayleighs 2 and we 
recall (Table [J) that W x {Ua) w 12.5 A in our corre- 
lation spectrum. Our result ([2"Tj) is close to the value 
/(Ha)/R = (0.129 ±0 015)(/ „(100 ^/MJysr" 1 ) found 
recently bv lWitt et all (|2010D . 

Averaging th e 100 ixm in t ensity from SFD and the Ha 
intensity from iFinkbeinerl (|2003[ ) over all regions with 



1R: 



10° 

4-7T 



photons cm 
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b > 60° gives a ratio 



/(Ha) 



R 



0.47 



J„(100/mi) \ 
MJysr- 1 / 



(22) 



The ratio of the coefficients in (|21|22[) gives a scattered 
fraction of Ha at high latitudes of 0.090/0.47 = 0.19. We 
add the statistical error of 5% in the Ha equivalent width 
and our conservative estimate of the uncertainty in the 
calibration shown in Figure |8] in quadrature, giving a 
scattered Ha fraction at high latitude of 



scattered Ha 
total Ha 



= 0.19 ±0.04 



(23) 



IWitt et all (|2010t) produced a scatter diagram of the scat- 
tered Ha fraction (their Fig. 6). The centroid of their 
distribution appears to be close to our measured value 
0. 19 ± 0.04 . 

iWood fc Reynolds! (|1999Q estimated that 5-20% of the 
Ha at high latitudes would be scattered light from H II 
regions. In a theo retical study of t he em ission spectrum 
of the diffuse Ha, IDong fc Draind (|2011[ ) concluded that 
the observed line ratios, including the low ratio of ra- 
dio free-free to Ha, could be understood if —20% of the 
diffuse Ha is actually reflected from dust, rather than 
emission from recombining gas in that direction. All of 
these results are consistent with our value 0.19 ± 0.04. 

We expect the scattered Ha to originate both from 
H II regions around young stars and from the WIM. Our 
line ratios seem to be more consistent with those of 
the W IM than with classical H II regions ( Madscn ct al. 
2006), which may indicate that the WIM dominates the 
line emission in our correlation spectra. However, this 
could also be due to the lack of nearby H n regions il- 
luminated by early O stars. Our line ratios do seem 
to be compatible with those in the 4> Per H II region, 
which is illuminate d by a B0.5+sdO system (Fig. 12 of 
IMadsen et aT][2006T) . 

5.4. Dust Luminescence 

We find no evidence of the Extended Red Emission 
(ERE) observed in some reflection nebulae (e.g., NGC 
7023) as a broad emission excess peaking near 7000A, 
with a FWHM -1500 A (jWitt fc Viihl IIM and ref- 
erences therein); a similar broad excess is reported 
to be present in the diffuse ISM ([Gordon et aliH998t 
iSzomoru fc Guhathakurtal H99fl IWitt et al.ll2008f) . The 
spectru m of light reflected by cirrus clo uds was mea- 
sured bv lSzomoru fe Guhathakurtal (|1998[) , who reported 
a broad peak in XI\ near 6500A that was attributed to 
luminescence, on the grounds that scattering is insuffi- 
cient. Our study also finds XI\ peaking near 6500A (see 
Figure [3]) but the observed spectrum appears to be con- 
sistent with what is expected for scattering. Indeed, our 
scattering models tend to cwer-predict the DGL around 
7000-9000 A (Figure O, depending on the adopted dust 
size distribution. 

If we posit that the slight excess emission between 5500 
A and 7500 A over our prediction with the ZDA04 dust 
model (see Fig. [3]) is an upper limit on the ERE, we infer 
that dust luminescence accounts for no more than 10% 
of the DGL in this wavelength range and less elsewhere. 
After calibrating the spectra and relating 7^(100 /im) to 



the total IR power using Eq. (fTTj) . this implies a ratio of 
ERE to infrared power of <0.4%, or a ratio of ERE to 
scattered power of <1%. 

The ERE will be a function of the intensity of illu- 
minating UV light and the dust column density. It is 
possible that our technique, which requires spatial vari- 
ations of the ERE correlated with 100 /im intensity over 
a —1° scale, is ill-suited to its detection. If the UV ISRF 
is particularly weak or spatially variable on small scales, 
it could at least partially explain the lack of ERE in our 
correlation spectrum. We have therefore performed two 
tests. 

First, we made our optical depth cutoff Ay < 0.5 more 
stringent to restrict the calculation to regions optically 
thin to UV photons. Reducing this cutoff by a factor of 2 
had a negligible effect on the correlation spectrum. Un- 
fortunately, when we reduced it further, the main stellar 
spectral features disappeared and new features similar to 
those in the high latitude spectrum of Figure|3]appeared. 
At such low levels of 100 /im emission, we expect extra- 
galactic contamination to be significant. 

We have also tried to estimate empirically how much of 
the small-scale variation (over ~1°) in 100 /im emission 
is due to variations in the illuminating starlight and how 
much is due to variations in the dust column density. 
We recalculated the correlation spectrum, correlating the 
SDSS intensities against dust optical depth rather than 
100 /im intensity. Our recovered correlation spectrum 
was nearly identical. Unfortunately, the resolution of 
the SFD temperature map used to convert from 100 /im 
intensity to optical depth, 0.7°, is too coarse to make this 
test conclusive. 

While we find no evidence of hidden ERE in our corre- 
lation spectrum, a more thorough investigation may re- 
quire instruments that can more directly probe the dust- 
scattered optical spectrum. 

6. SUMMARY AND CONCLUSIONS 

In this paper, we have measured the spectrum of the 
scattered component of the DGL using 90,000 calibration 
spectra from SDSS by correlating against 100 /im inten- 
sity measured by COBE and IRAS and reduced by SFD. 
The correlation spectrum of the DGL is consistent with 
scattering of photons emitted by stars and by ionized 
gas. Its continuum and lines show interesting variations 
across the sky, which could be due to differences in the il- 
luminating starlight, in the dust properties, or structure 
in the Galaxy. Our spectrum is not calibrated because 
the correlation neglects (unknown) errors and small-scale 
structure in the 100 /jm intensity. However, an asymp- 
totically unbiased estimator agrees well with simplified 
radiative transfer calculations and allows us to calibrate 
the average spectrum of the DGL. In calibrating the data, 
we also provide an estimate of the errors and small-scale 
structure in the SFD 100 /im map. 

With a simplified radiative transfer calculation assum- 
ing a plane-parallel galaxy, we show that the spectrum 
of the DGL can discriminate between dust models. The 
ZDA04 dust model is preferred to the WD01 model, 
which has more large grains and scatters too much red 
light. The ZDA04 model fits the data very well without 
requiring dust luminescence as a source of ERE. Indeed, 
our radiative transfer calculations tend to over-predict 
scattering redward of 7000 A. 
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Our measurement of Ha in the DGL allows us to in- 
directly measure the fraction of scattered Ha at high 
Galactic latitude. We find an avera ge value of 0.19 ± 
0.04, co nsistent wi t h the results of iWood k, Reynolds! 
(11999ft. IWitt et all (12010ft and the value inferred by 
iDong & Draind (12011ft . The line emission also allows us 
to constrain the properties of local H n regions. The lack 
of emission from species with high ionization energies is 
consistent with the lack of nearby early O stars, and 
the average scattered spectrum is consistent with gas at 
-7500 K. 

Our results are the product of the extraordinary SDSS, 
the only dataset of its kind. It would be possible, though 
difficult, to incorporate the new Bary on Oscillations 
Spect roscopic Survey spectra (BOSS, lEisenstein et all 
12011ft into the analysis. The new fibers are smaller (2" 
instead of 3"), losing far more light and making spec- 
trophotometric calibration much more difficult. The 
DGL is a sufficiently weak signal that it requires thou- 
sands of blank sky spectra to extract useful results. 
These spectra also need to be in regions with significant 
dust columns, regions that tend to be avoided by large 
extragalactic surveys. The SDSS spectra will likely re- 
main the best tool for analyzing the DGL for some time. 

While the optical data are unlikely to improve in the 
near future, WISE and AKARI may allow significant im- 
provements to the infrared map. Our results indicate 



that the measurement errors and small-scale structure 
in the SFD 100 jxm map are nearly as large as the ~1° 
variations to which our correlations are sensitive. By 
combining high resolution infrared data with the SFD 
maps, we could improve the quality of our DGL spectra 
and resulting constraints on Galactic stars and dust. 
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APPENDIX 

SCATTERING AND ABSORPTION IN A PLANE-PARALLEL EXPONENTIAL GALAXY 

We test dust models by calculating scattered spectra in an infinite plane-parallel galaxy, using vertical distributions 
of stars and dust appropriate to the Solar neighborhood. The Sun is assumed to lie in the midplane, z — 0. We 
calculate the absorption (and hence, infrared emission) and scattering along each line-of-sight, neglecting multiple 
scatterings. 

We denote the (wavelength-dependent) optical depth of dust above height z along a vertical line-of-sight as t\{z), 



a exi (X)p(z') dz' 



(Al) 



where a cx t is the extinction cross-section and p is the dust density. Consider a sheet of stars with uniform surface 
power density Ti\ at height z s . The optical depth from an annulus of stars at a distance R from the grain's projected 
position in the stellar sheet to a grain at height z is 



A\(z,z 3 ,R) = \t x (z) -t\(z s 



v /( z _ Zs )2 +jR 2 



(A2) 



The total flux density (neglecting multiply scattered photons) incident on a grain at height z is then 

Fx(z, Zs ) = j Q 2*RdR 47T[{z _ Zs)2+R2] 

E A^ Ei (\t x (z)-t x (z s )\) , (A3) 

where E\ is the exponential integral and t\ is re lated to z by Equation (|Al|) . The total reradiated intensity is an 
integral of the absorbed fraction of Equation (|A3[) over all wavelengths and stellar sheets. Defining ui\ as the albedo, 
the total infrared intensity Jtir from a sightline at Galactic latitude b is 



esc |6| 
8tt 



(1 - ux) dX 



dz s 



■ dz* 



rx(0) 



dr' E x (\t'-t x (z s 



(A4) 



A dust model, like the iDraine &: Lil (|2007ft model used in Equation (IT71) . is needed to convert Jtir into a specific 
intensity I u at infrared frequency v. 

Anisotropic scattering makes it more difficult to calculate the optical intensity. Consider an annulus of the stellar 
sheet centered directly below a dust grain, and let 9 sweep out the annulus, with 9 = pointed away from us. If the 
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dust grain is along a sightline at Galactic latitude 6, the law of cosines gives the required scattering angle £ as 

(z — z s ) 2 — Rz cos 8 cot b , . „. 

cos£ = v ; . (A5) 

^(z 2 cot 2 b + (z - z s ) 2 ) (R 2 + (z- z s ) 2 ) 
The intensity of light scattered in our direction from a stellar sheet at height z s is then given by 

f°° f 27V , s s a exp [- A x (z,z s ,R)\ , 

h, sca (z,z s ,b)= RdR dO^x cosO . * ' °' J , A6 

Jo Jo 47r((z-z s ) 2 + i? 2 ) 



with 0a being the normalized phase function and A\ dehned by Equation (|A2j) . We use a H enyey-Greenstein phase 
function for simplicity. If <j> is constant, the 6 integral is trivial and we recover Equation (|A3[) modulo a factor of 4-7T. 
Scattered light is further attenuated by dust along the line-of-sight. Neglecting multiple scatterings, we finally have 

^a. sca(2 s , o) = wacsc |6| / exp [- esc | b\ (r A (0) - r )J dr / RdR d6> </>a (cosf) — ,, . 2 2 . - , 

JO JO JO 47r U z z sj + A ) 

(A7) 

where z a nd ta are related by Equation (jAll) and A a is defined by Equation (| A2|) . As for Equation (IA4[) , we integrate 
Equation (|A7I) over z s . We do not integrate over wavelength; the spectrum of Ja , sca is shown in Figure [3] normalized 
by the infrared intensity vl v at 100 /zm, which is related to Itir by Equation |T7|) . The ratio /A.sca/^TiR depends 
very weakly on Galactic latitude b, the normalization of the optical depth, and the details of the stellar and dust 
distributions. Once an ISRF and dust model are supplied, there are no other free parameters. 
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